Scale invariant properties of public debt growth 
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Public debt is one of the important economic variables that quantitatively describes a nation's 
economy. Because bankruptcy is a risk faced even by institutions as large as governments (e.g. 
Iceland), national debt should be strictly controlled with respect to national wealth. Also, the 
problem of eliminating extreme poverty in the world is closely connected to the study of extremely 
poor debtor nations. We analyze the time evolution of national public debt and find "convergence" : 
initially less-indebted countries increase their debt more quickly than initially more-indebted coun- 
tries. We also analyze the public debt-to-GDP ratio 1Z, a proxy for default risk, and approximate 
the probability density function P(1Z) with a Gamma distribution, which can be used to establish 
thresholds for sustainable debt. We also observe "convergence" in 1Z: countries with initially small 
7Z increase their 1Z more quickly than countries with initially large 1Z. The scaling relationships 
for debt and 7Z have practical applications, e.g. the Maastricht Treaty requires members of the 
European Monetary Union to maintain 1Z < 0.6. 
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I. INTRODUCTION 

Just as an individual is expected to control his/her 
debt to asset ratio, so is a government expected to control 
its national debt as a function of the country's wealth, 
measured e.g. by its gross domestic product (GDP). In 
a dynamic global economy, excessive borrowing cannot 
persist indefinitely, as creditors are bound to call in large 
loans. While a country suffers financial problems when 
its GDP does not increase fast enough, even more seri- 
ous trouble begins when its debt increases faster than its 
GDP. While national GDP has been the topic of many 
studies on economic growth [U-IH, the empirical analysis 
of public debt has lagged due to lack of comprehensive 
data. 

Large sets of public debt data, dating back several 
decades, and ranging from poor to rich countries, have 
recently become available. Here we use concepts of statis- 
tical physics to analyze public debt data for a wide cross- 
section of economies including underdeveloped, develop- 
ing, and developed countries. The total public debt data, 
along with total GDP data, are available at the Inter- 
American Development Bank and are compiled and 
analyzed in Refs. @, 0|- Population data are available 
by the World Bank, and can be reconstructed through 
GDP and per capita GDP data compiled in Ref. [8| . We 
deflate all USD amounts to in units of the USD in the 
year 2000 @. In our analysis, we compare public debt 
only within the same country, in order to avoid any dif- 
ferences in the theoretical and practical definition of debt 
and the reporting of debt by various countries, an issue 
pointed out in Refs. @, 0)E3l- O ur results are robust 
with respect to mis-reporting and ambiguous definitions 
of public debt [15. In Fig. [Owe plot the debt-to-GDP 
ratio 7Z for many countries, grouped in panels (a-d) by 



common historical, geographical, and financial factors. 
In Fig. QJd) we plot the average debt-to-GDP ratio for 
three subgroups corresponding to World Bank Income 
Group (IG) classifications, and observe relatively high 
levels of TZ among the poorest countries. 

Economic growth theories predict that GDP should 
"converge" towards equality, with wealthy countries ex- 
periencing smaller relative growth rates than poor coun- 
tries. However, the opposite has been found for economic 
wealth data jllTfl3l ]. So we address the question, what 
are the growth dynamics for public debt? To answer this 
question, we analyze a comprehensive database of na- 
tional public debt and GDP to investigate the dynamics 
of debt growth and growth in TZ. 

With the current global credit crunch, and several no- 
table recent national defaults, it is important to address 
sustainable public debt, defined as the amount of debt 
where the receiving country is capable of meeting its 
current and future debt obligations [3, [l5[ . The total 
current government debt D(t) increases from last year's 
debt D(t — 1) partially due to interest payments on the 
debt D(t — 1) at interest rate Id if), and partially because 
of the current primary deficit, defined as the difference 
between spending S(t+ 1) and taxes T(t + 1) [16]. Thus 

D(t) = [l+I D (t- l)]D(t - 1) + [S(t) - T(t)} . (1) 

We consider three possible scenarios for public debt 
growth dynamics: 

(i) Growth rates of the country debt do not depend on 
the initial debt level. 

(ii) A more indebted country has a larger debt growth 
rate than a less indebted country, so that relative 
differences between debt across countries increases 
over time (divergence). 
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(iii) A more indebted country has a smaller debt growth 
rate than a less indebted country, so that relative 
differences between debt across countries decreases 
over time (convergence) . 

These three scenarios have different implications for in- 
vestors, who will only accept government debt up to some 
ceiling. Hence, one would expect that more indebted 
countries would increase their debt more slowly than less 
indebted countries. 



II. EMPIRICAL RESULTS 

To ascertain which of the three debt scenarios is better 
supported by empirical facts, we define for country i the 
annualized logarithmic growth rate of per capita initial 
debt di(t) between years t and t + At 

n(t,t + At)^ l -2iM±^m. (2) 

We compare ri(t, t + At) to di(t), assuming ri(t, t + At) 
depends on debt size by 

n(t,t + At) *±a-p \og[(k{t)] . (3) 

The functional form of Eq. ((3]) can also be expressed as 

\og[di(t + At)} = a At + (1 - 0At) log[di(t)] (4) 

If ft > 0, there is convergence in per capita debt data 
across countries, since initially more indebted economies 
tend to increase their debt slower (smaller r*j(£, t + At)) 
than initially less indebted economies. Hence, (3 repre- 
sents the "speed of convergence" , a concept introduced 
for per capita GDP data in Ref. [HJ]. A larger posi- 
tive value of f3 results in faster convergence, equalizing 
the per capita debt across all countries more quickly. If 
j3 < there is divergence in debt data, where initially 
less- indebted countries with smaller di(t) increase their 
debt slower than initially more-indebted countries. 

In Fig. [5] we plot the 1990 per capita debt of more 
than 80 countries, representing low-, medium-, and high- 
wealth countries, considering several relationships. First, 
we compare the per capita debt over (a) 15-year and (d) 
7- year time horizons. We find the slope S = (1 — (3 At) of 
the regression in Eq. (Q} is less than one, requiring (3 > 
which corresponds to scenario (iii). 

To confirm the convergence across countries for other 
time horizons, Fig. [3] shows the value of S for varying 
initial d(t) and time horizon At in Eq.Q. We find S = 
1 — (3 At < 1 for most horizons At, implying convergence, 
where less-indebted countries increase their debt faster 
than more-indebted countries. However, there is a period 
in the beginning of the 1990's that is the exception, with 
S > 1 and /3 < 0. This period of divergence in per capita 
debt may be related to the 20-year lows in interest rates 
which may have resulted in increased borrowing, even 
among heavily indebted countries. Since 1995, the values 



of S have returned to values less than one, indicating a 
return to convergence. 

We also analyze the growth rates of per capita GDP 
and confirm the divergence across countries observed 
originally in [TlT[l3j . In Fig. [4] we plot for per capita 
GDP, the analogous regression S values that we plot in 
Fig. [3] for per capita public debt. For all periods At and 
initial years analyzed, we find values of S > 1 indicating 
divergence. 

A natural question is - how does the per capita debt 
di vary across all countries and by income group? Power 
law probability density functions (pdf) have been ob- 
served for total country GDP [17j and per capita GDP 
fl8j |. Fig. [5] shows the pdf P(d) for all countries an- 
alyzed over the 36-year period 1970-2005. We observe 
large variations across income groups, where low income 
countries typically have relatively small per capita debt 
values reflecting their small per capita borrowing capac- 
ity. In contrast to the zipf-rank curves for GDP with 
Cgdp ~ 1 corresponding to pdf scaling exponent expo- 
nent 1 + 1 /Cgdp ~ 2[17|], we observe in Fig. [Sunset) a 
scaling value Cd ~ 0.3 corresponding to a relatively large 
pdf scaling exponent 1 + l/Cd ~ 4.3. 

In a country where both GDP and debt grow with 
time, one must analyze the dynamics of both debt and 
GDP. Since a debt that is large for Luxembourg is not 
large for the U.S., various indices have been proposed 
in order to compare the burden of debt to the ability of 
the country's economy to generate income. These include 
1Z pjj [2CJ, so we apply the convergence analysis of Eq. (TJ| 
to IZ(t) obtaining 

log[R 2 (t + At)} = a' At + (1 - /3' At) log[7^(i)] . (5) 

Fig. [2] compares lZ(t) over (b) 15-year and (e) 7-year 
time horizons. Fig. [7] shows S' = (1 — f3'At) < 1, im- 
plying convergence f3' > 0, over a large range of At-year 
horizons for initial year t. 

A responsible government is expected to monitor si- 
multaneously the growth of debt and GDP [2l| . By bor- 
rowing money, a country may increase TZ(t) for some 
time, but clearly lZ(t) cannot increase indefinitely, as 
increased debt can negatively affect GDP growth |22j . 
Banks prefer individuals with large incomes and small 
debts. Banks also prefer countries that have, for a given 
GDP level, small relative debt. Fig. [TJ provides the an- 
nual trend of lZ(t) for several groups of countries with 
common geo-politial backgrounds. 

Debtor default risk is estimated by many rating agen- 
cies and financial organizations. 7Z(t) is an important 
quantity for determining the ability of a debtor to make 
debt payments. For large lZ(t) there is a larger proba- 
bility that the debtor will not be able to make timely 
payments or be able to prevent further debt increase 
with time, scenarios that lead to credit default. In or- 
der to quantify the risky debt levels, we collect the TZ(t) 
values of all countries analyzed over the 36-year period 
1970-2005 and plot the pdf P(K) in Fig. E We find 
< 1Z >= 0.57 ± 0.54, and we fit the pdf to a Gamma dis- 
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tribution P(TZ) oc exp[-TZ/TZ c ] with k = 2.0 ± 0.1 

and 1Z C = 0.30 ± 0.01, using the maximum likelihood es- 
timator. The extreme value statistics of Gamma pdf can 
be used to define thresholds for sustainable debt. 

In order to analyze the countries that have large lZ(t) 
and a high risk of default, the countries which consti- 
tute the pdf tail, we plot rank-frequency curves in Fig. 
|5J inset). The Zipf plots show a power law over three 
orders of magnitude, with scaling exponent Ctz ~ 0.4 
corresponding to P(1Z) ~ 1Z~ 3 5 . 



III. MODEL 

Our analysis, performed across a wide cross section of 
countries, confirms the existence of convergence in public 
debt. This is opposite of what is in GDP data [lll - ll3j . 
where the speed of convergence /3 is negative. We now 
discuss how to model the scaling result we obtain, and 
how to use the scaling result obtained for GDP and public 
debt. Figs. [2{c) and [2jf) compare the per capita debt 
to per capita GDP for the years 1990 and 2005. The 
typical relationship between debt and GDP shows a scale 
invariant form, 



g~A<T 



(6) 



where g is the per capita GDP and d is the per capita 
debt. In Fig. [7] we plot the values of 7 for the set of 
countries analyzed in each yearly data set. 

In order to model debt dynamics, we assume that the 
functional dependence in Eq. (j6|) is time invariant. Note 
that if g(t) and d(t) grow exponentially with different 
growth rates, r g and r^, the relationship between g(t) 
and d(t) still has the form of a power law, with r g = 7 r d . 

We may consider the dynamics of public debt by as- 
suming that the government borrows B(t), a fixed pro- 
portion of GDP given by B(t) = D(t) - D(t - 1]_= 
AD(t) — b G(t), with deficit ratio b = constant |23j . 
Thenrr, = AD/D = b G(t)/D(t), and 



r d = b g(t)/d(t) - r p 



bN/d{t) 



1-7 



(7) 



where r pop denotes the population growth rate [2C 
Hence, 



dr d _ -bN(l - 7) 



dd 



d 2 -> 



l/d 2 -~< 



(8) 



We observe from Fig. [7] that 7 < 1 so that dr^/dd 
— l/d. For this reason, we use the regression r c i = a 



(3 log[df(t)], which agrees well with the data in Figs. ([2]) 
and jS|). 



IV. DISCUSSION 

In summary, we demonstrate convergence in per capita 
public debt across a wide set of countries during the pe- 
riod 1970-2005, a result of general interest for complex 
systems researchers as well as for creditors. Our anal- 
ysis is made possible by new comprehensive data sets, 
and extends empirical surveys previously performed on 
country GDP which found divergence in country GDP. 
While divergence in country GDP implies that economic 
wealth is moving away from global equality, convergence 
in per capita debt implies that indebtedness is becoming 
an economic standard. Furthermore, convergence in the 
1Z implies that relative differences in indebtedness across 
countries is also decreasing over time. Some economists 
believe that convergence across all countries is possible 
through globalization fPH, [24T - [26| and access to open mar- 
kets [27| ■ While public debt can be used to invest in a 
country's development via physical infrastructure, tech- 
nology, and social programs, its use requires responsible 
governance. Corruption [H, [2!| and the misuse of pub- 
lic debt can lead to insurmountable debt contributing 
to financial crisis, which can cause further increase in 
debt levels through exchange rate depreciation (3Tj, [3l| . 
There are also instances of extremely poor debtor na- 
tions that are unable to meet their current and future 
debt obligations (See Fig. [T|). Recent programs such as 
the Heavily Indebted Poor Countries (HIPC) Initiative, 
sponsored by the World Bank and the IMF, and the Ju- 
bilee 2000 Campaign-to-Drop-the-Dept, have called on 
debt cancellation for extremely poor debtor nations as a 
crucial step in the UN Millennium Project to eliminate 
extreme poverty (32Tj34l |. Further, debt has become a 
problem for not only the extremely poor countries. With 
the current global credit crunch, their is an increased 
need for responsible use of government debt. 
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Year Year 

FIG. 1: Illustration of debt-to-GDP ratio 7Z(t) for (a-d) several countries and (e) the average trend for three subsets according 
to World Bank Income Group (IG) classifications, (a) Countries with turmoil, as in the case of Israel and Sudan which 
experienced periods of war, and Guyana which formed a new government, can experience periods of extremely high 1Z. (b) 
The Argentine financial crisis of the late 1990's placed many South American countries in a state of debt stress, (c) The debt 
burden in many African countries, some which are also plagued by war and political turmoil, is increasing with time. Many 
of these countries are candidates for foreign debt cancellation (jubilee), (d) Several large countries, (e) The average 1Z is 
increasing with time, with low-income countries currently at relatively high levels reflecting the burden of unsustainable debt, 
(f) Mediterranian countries with recent concern of default risk. 
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FIG. 2: The 3 possible relations between per capita debt, per capita GDP, and 1Z for a set of low-income, medium-income, 
and high-income countries. Regression of per capita debt d(t) in 1990 versus per capita debt d(t + At) for (a) At = 15 years 
and (d) At — 7 years later indicates that there is convergence in the per capita debt data across the countries over the larger 
time horizon in panel (a) with S = 1 — /3At < 1 implying f3 > 0. Regression of lZ(t) in 1990 versus 1Z(t + At) for (b) At — 15 
years and (e) At = 7 years indicates that there is also convergence over these time horizons. (c,f) We also illustrate the scale 
invariant relation between between per capita debt and per capita GDP in units of 10 3 USD per person in the year 2000. 
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FIG. 3: The slope S = 1 — /3At of the log-log regression in Eq. Q is typically less than one for a given time horizon of At 
years after the initial year t labeling each panel. Values of S < 1 imply convergence (f3 > 0) in per capita debt over the Ai-year 
period. Ft 2 values for log-log regressions are typically > 0.7 for 10 < At < 20 years, and are typically > 0.8 for 1 < At < 10 
years, while the number of countries for each regression varies between 22-25 (1970) and 103-120 (1995). Interestingly, we 
observe a short period of divergence, where for initial years 1980 and 1985, the value of S > 1 for At ~ 12 years corresponding 
to the mid-1990's. 
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FIG. 4: For comparison, we also analyze GDP growth rates and find divergence for all periods analyzed. The slope S of the 
log-log regression in Eq. Q for per capita GDP is typically greater than one for a given time horizon of At years after the 
initial year, which we label in each panel. Values of S > 1 imply divergence (/3 < 0) in per capita GDP over the At-year period, 
consistent with the seminal work of fH. [T^|. 




per capita debt , d 

FIG. 5: Right-skewed probability density function P(d) of per capita debt d for all countries analyzed over the 36-year period 
1970-2005 in units of 10 3 USD per person in the year 2000. (Inset) Zipf rank-frequency plot for d values for all countries and for 
three subgroups according to World Bank Income Group (IG) classification. Countries in the high and medium IG categories 
contribute to the wide range of d values, whereas poor countries have relatively small d values reflecting their small per capita 
borrowing capacity. We also show a power-law with exponent £ = 0.3 for comparison. 
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FIG. 6: The slope S' = 1 — /3' At of the log-log regression in Eq. ((5]) of 1Z. We find S' < 1 a given time horizon of At years after 
the initial year t labeling each panel. Values of S < 1 imply convergence (/3 > 0) in per capita debt over the At-year period. 
We plot values of S' only for log-log regressions with R 2 > 0.3. The number of countries for each regression varies between 
22-25 (1970) and 103-120 (1995). 




FIG. 7: Annual trend of power-law exponent 7, which quantifies the scaling relationship between debt and GDP as g ~ A <P 
in Eq. (0. We find 7 < 1 which reflects the burden of debt on country GDP. 
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R = D/GDP 

FIG. 8: Probability density function P(1Z) for all countries analyzed over the 36-year period 1970-2005. For comparison, we 
approximate P{7V) with a Gamma distribution with parameters k — 2.0 and 1Z C = 0.30. The extreme statistics of P(TZ) can 
quantify the threshold of sustainable debt. (Inset) Zipf rank-frequency plot for 1Z values. We separate data into three subgroups 
by World Bank Income Group (IG) classification to demonstrate the range of large 1Z values within different wealth groups. 
The straight line is a power law with exponent £ = 0.4 for comparison. 



